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We calculate exactly the quantum mechanical, temporal second-order coherence function for a 
single-mode, degenerate parametric amplifier for a system in the Gaussian state, viz., a displaced- 
squeezed thermal state. The calculation involves first the generation of the Gaussian state and 
subsequent measurements of two photons a time r > 0 apart. The generation of the Gaussian state 
by the parametric amplifier insures that the temporal second-order coherence function depends only 
on r. 

PACS numbers: 42.50.-p, 42.50.Ct, 42.50.Ar, 42.50.Xa 


Quantum optical systems give rise to interesting nonclassical behavior such as photon antibunching and sub- 
Poissonian photon statistics owing to the discreetness or photon nature of the radiation held jlj. These nonclassical 
features can best be quantified with the aid of the temporal second-order quantum mechanical correlation function 
g( 2 \r) and experimentally studied using a Hanbury Brown-Twiss intensity interferometer modified for homodyne 
detection )2j. 

Consider the unitary transformations 

S(—t;)D{—a)aD{a)S{£,) = a cosh r — a) exp (id) sinhr + a, (1) 

and 

S{—t;)D{—a)a^D{a)S{t;) = a) coshr — aexp(— iO) sinhr + a*, (2) 

with the displacement D{a) = exp (cred — a*a) and the squeezing S(£) = exp ( — | a) 2 + 4ya 2 ) operators, where a 
(cr) is the photon annihilation (creation) operator, £ = r exp {id), and a = |a| exp (iip). 

The operators {aa^ — a*a) and (—|ad 2 + 4j-a 2 ) do not commute with their commutator and so the product 
S{—£)D{—a) cannot be expressed in the exponential form exp{iHt/h). However, we can find an operator H such that 
the Eqs. (1) and (2) are reproduced. The Hamiltonian for degenerate parametric amplification, in the interaction 
picture, is 

H = ca) 2 + c*a 2 + ba + b*a^. (3) 


One has that 


exp(ziTr/?i)aexp(— iHr/Ti) = a cosh(2|c|r/?i) — ia* exp {i\) sinh(2|c|r/?i) — ^* s ^ n ^ 1 ^‘ J l c l r /^) 

2|c| 


+ b exp (ix) 


cosh(2|c|r/?i) — 1 

2M 


where c = |c|exp(zy). Comparing Eqs. (1) and (4) yields 


£ 0 ) 



r(r) exp {id), 


where 


r(r) = 2\c\t/H, 


( 4 ) 

( 5 ) 

( 6 ) 
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and 
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a[r) = —ib * 


sinhr(r) 

2H 


ibexp (iO) 


coshr(r) — 1 

2|c| 


(7) 


Note that r(0) = a(0) = 0. 

The quantum mechanical temporal second-order coherence [lj is given by 

(2 ) = (aHt)a\t + T)a(t + T)a(t)) 

(ai(t)d(t))(at(t + r)a(f + r)} ’ 

with t > 0, where the operators are in the interaction picture and the expectation values, representing an average 
over the initial photon states, are determined by the thermal density matrix 

Po = exp (—/37iu;n)/Tr[exp (— /3hion)], (9) 

where h = a)a, /3 = 1 /(fcsT), and n = Tr[/3oro]. In quantum optics, one usually deals with statistically stationary 
fields, viz., the correlation functions of the field are invariant under time translation and so the correlation functions, 
for instance, the temporal second-order coherence (8), would be a function of the single variable r. Note, however, this 
requires that the thermal density matrix (9) commutes with the Hamiltonian (3) that governs the time development of 
the system. This, however, is not so. Nonetheless, we have suppressed the dependence on t in the two-time correlation 
(8) since g^ 2 \r) can be expressed as 


( 2 ) / \ = (d t (0)d t (r)d(r)d(0)) 
(at(0)a(0))(at(r)a(r))’ 


( 10 ) 


where all the expectation values in (10) and henceforth are traces with the Gaussian density operator, viz., a displaced- 
squeezed thermal state, given by 

p G = D(a)S(Z)p 0 S(-Z)D(-a), (11) 


where a = a(t) and £ = £(t). Accordingly, the system is initially in the thermal state po. After time t , the system 
evolves to the Gaussian state pc and a photon is annihilated at time t, the system then develops in time and after a 
time r another photon is annihilated. Therefore, two photon are annihilated in a time separation r when the system 
is in the Gaussian density state pc- 
Now 


(a t (r)a(r)) = {n + 1/2) cosh[2(r + r(r))] - 1/2 + |A(r)| 2 , (12) 

where 

A{t) =acoshr(r) — ia* exp (i0) sinhr(r) + a(r), (13) 

with r(r) and a(r) defined by Eqs. (6) and (7), respectively. One has for r = 0 that 

(al(0)a(0)) = (n + 1/2) cosh(2r) — 1/2 + |a| 2 . (14) 

The calculation of the temporal second-order correlation function (10) gives 

( 2 ), \ _ 1 rt 2 (r) + s 2 (r) + [aA* (t) + a* A(t)]u(t) - [aA(r) exp (-id) + q*A*(r) exp (i0)] s(r) 

9 T (at(0)a(0))(at(r)a(r)) 

where 

n(r) = (n + 1/2) cosh (2 r + r(r)) — (1/2) coshr(r), (16) 

and 

s ( t ) = {n + 1/2) sinh (2 r + r(r)) — (1/2) sinhr(r). (17) 

Result (15) reduces to the known results for r = 0, for the vacuum state with h = 0 @ and for the thermal state with 
n ^ 0 3, since r(0) = 0 and A(0) = a. 
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The parameters c and b in the degenerate parametric Hamiltonian (3) are determined by the parameters a and £ 
of the Gaussian density of state (11) via 


tc = —i—r exp (id) 


(18) 


and 


tb = 



^aexp (—id) + a* coth(r/2)^r, 


(19) 


where t is the time that it takes for the system governed by the Hamiltonian (3) to generate the Gaussian density of 
state pc from the initial density of state po- 

Note that the ratio b/\c\ from (18) and (19) depends only on the parameter a of the displacement operator D(a) 
and the parameter £ of the squeezing operator S(£). Hence, a(r) defined by (7) depends on r only through r(r), 
which is defined by (6). Therefore, A(t) defined by (13) depends on the parameters a and £ and on the variable r 
only through r(r) also. Accordingly, <p 2 )(r) depends on the parameters a and £ of the density pc and on the variable 
r only through r(r). Therefore, for given a and £, the temporal second-order coherence function g^(r) depends on 
r only on r(r) and the mean photon number n of the thermal state po . 
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